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Abstract 

We perform the scanning tunneling spectroscopy based superconductor-vacuum-superconductor ana- 
logue to the seminal McMillan and Rowell superconductor-insulator-superconductor device study of 
phonons in the archetypal elemental superconductor Pb [W. L. McMillan and J. M. Rowell, Phys. Rev. 
Lett. 14, 108 (1965)]. We invert this spectroscopic data utilizing strong-coupling Eliashberg theory to 
obtain a local a 2 F(co) and find broad underlying agreement with the pioneering results, highlighted by 
previously unobserved electron-hole asymmetries and new fine structure which we discuss in terms of both 
conventional and unconventional superconducting bosonics. 

PACS numbers: 74.55.+V, 74.25.Kc, 74.45,+c, 74.78.-w 
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It has been well established that the Cooper pairs of conventional superconductors are bound 
by a "pairing glue" mediated by phonons 03I21- The verification of phonon-mediated pairing was 
facilitated by a combination of experimental evidence including an isotope effect study on ele- 
mental Hg 0, tunneling studies on Pb devices by McMillan and Rowell (M&R) and Giaever et 
al. flU SI, and phonon measurements on Pb by Brockhouse et al. [5J. Since tunneling studies 
played such a crucial role in the original identification of the pairing glue, point contact and de- 
vice tunneling methods have been in continual employment to study the bosonic structure of both 
conventional and unconventional superconductors alike |0-[8l. Additionally, with recent techno- 
logical improvements in scanning tunneling spectroscopy (STS) experiments, bosonic informa- 
tion has been recorded for various high transition temperature (high-r c ) superconducting systems 
utilizing metallic tips [1914151. This STS-based setup proves superior to superconductor-insulator- 
superconductor (S-I-S) device tunneling for its ability to measure spectra at single atomic locations 
rather than large spatial averages, in addition to STS utilizing the cleaner vacuum barrier over the 
device's sandwiched insulator. 

Remarkably, despite the increasing prevalence of STS bosonic studies, no control experiment 
has ever been completed where vacuum-based STS is performed on a known conventional sys- 
tem to measure and fully invert the bosonic substructure, verifying it against the known result. 
This is especially important since certain details of the tunneling process are still under inves- 
tigation [fT6l - l22l . Though conventional superconductors are better understood than their uncon- 
ventional counterparts, acquiring high-resolution STS spectra on conventional superconductors is 
more challenging. This difficulty stems from significantly smaller gap sizes (an order of magni- 
tude) which lead to more stringent requirements on noise minimization for resolution comparable 
to that of high-r c systems. (These requirements become even more stringent when recording a 
higher harmonic of the modulation voltage to acquire d 2 I/dV 2 ). There have been numerous prior 
results for vacuum STS on conventional superconductors but these results have either ignored the 
bosonic substructure of the spectra while focusing on in-gap states, or not possessed the required 
energy resolution and noise levels for clean spectroscopy [|23l - l25l . In light of this, as diagrammed 
in Fig. 1(a) we studied STS on the archetypal elemental superconductor Pb, employing a super- 
conducting Pb tip for feature sharpness and to construct a true vacuum analogue (S-Vac-S) to the 
M&R device study (S-I-S) JTJ. (The reasons for this specific material choice, as well as the de- 
tailed experimental method, can be found in the Supplemental Material [SM]). We believe this 
to be the first inversion study of the phonon spectrum of a conventional superconductor utilizing 
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FIG. 1: (a) A diagram of the experimental S-Vac-S setup, (b) Raw dl/dV spectrum at a single location. The 
negative sample biases, corresponding to occupied electronic states, are colored in blue while the positive 
biases, corresponding to unoccupied electronic states, are colored in magenta. The spectrum in gray is the 
normal state spectrum acquired in a magnetic field B = 1 T. (c) d 2 I/dV 2 by numerical differentiation of 
the data in (b). The yellow band highlights the range of interest which we zoom in on in (d), showing 
high-resolution d 2 I/dV 2 acquired by HHS. 
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ultra-high vacuum STS in true tunneling mode. 

To verify the superconducting nature of the sample and tip we measured a dl/dV spectrum 
using the standard lock- in technique [Fig. |T|[b)]. The spectrum shows the familiar structure for 
tunneling between two inequivalent superconductors with tall coherence peaks located at A = 
A t ip + A sample (A sample > A tip ) and two in-gap features at A' = A sample - A tip . From this we obtain 
A = 2.0 ± 0.02 mV and A' = 0.42 ± 0.02 mV which corresponds to superconducting energy gaps 
Asample = 1.21 ±0.02 mV and A t i p = 0.79 ±0.02 mV. The A samp i e value agrees well with the 
expected bulk value and the A t j p value is consistent with other experimental gap measurements on 
few-monolayer Pb systems Il26ll27ll . 

Continuing outwards from the coherence peaks, sharp dips are observed in dl/dV on both 
sides of the Fermi energy £p (V = 0) at energy £dip = 2.72 ± 0.02 mV. This value is equal (within 
experimental error) to A samp i e ± 2A t i p , but this may be a coincidence. Following the large dips 
are several step like features that are well known to be associated with electron-boson coupling 
(phonons, in this case). These features are best analyzed by examining d 2 I/dV 2 vs. V [Fig. 1(c)] 
which can be directly related Q~| to the effective electron-boson coupling function a 2 (co) and the 
boson density of states F(co) for bosons of energy ft). We will discuss the dl/dV dips in detail 
later in the Letter and first focus on the bosonic features (which are reported in terms of Q. where 
Q = ft) -A). 

The largest d 2 I/dV 2 signal [Fig. 1(c)] comes from the sharpness of the coherence peak which 
is two orders of magnitude larger than the other peaks. In order to resolve the much weaker signal 
from the boson substructure it is necessary to tune the settings of the tunnel junction and lock-in 
preamplifier to zoom in specifically on these higher energy features. To highlight these complex 
spectral features we utilize a second lock-in preamplifier tuned to twice the modulation frequency 
to directly resolve d 2 I/dV 2 , which we will refer to as higher-harmonic spectroscopy (HHS) [Fig. 
1(d)]. 

In order to compare this data with the M&R device data, the M&R result has been digitized, 
thermally broadened, and overlaid with this study's result (see [SM]). It can be seen that the 
STS based S-Vac-S data is in excellent agreement with the M&R device S-I-S study (Fig. [2]). The 
differences between the two sets of data will be highlighted and further analyzed later in the Letter. 

Unlike HHS, the M&R method calls for taking dl/dV in the superconducting state, dividing it 
by dl/dV in the normal state, and differentiating the result. This was necessary in order to divide 
out any spurious signals resulting from the insulating tunneling barrier which might contribute 
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FIG. 2: Smoothed d 2 I/dV 2 vs. D. comparing the occupied and unoccupied spectra from this experiment 
with the M&R device result (H. (Occupied curve is —d 2 I/dV 2 ). The M&R data has been thermally 
broadened from the 0.8 K at which the data was recorded to 4.2 K to match this experiment (see [SM]). The 
occupied electronic states spectrum is accompanied by the raw data (black circles). 

inelastically to the tunneling process. Since in our experimental setup the vacuum acts as a perfect 
insulating barrier, both methods should give similar results. In order to substantiate this claim 
we repeat this experiment using the traditional M&R method utilizing a 1-T magnetic field (the 
critical field of Pb is ps 0.080 T Il28"1 ) and compare the results. It can be seen that the normal 
state curve is essentially flat when compared to the superconducting state [Fig. [TJb)], so it can be 
accurately surmised that dividing dl/dV |sc by an approximately fiat line will result in essentially 
the original curve. For thoroughness we nevertheless complete the M&R-prescribed recipe to 
obtain d 2 I/dV 2 using the superconducting and normal state data of Fig. 1(b). The resultant spectra 
is overlaid on the HHS result in Fig. ^a). It is evident that the spectra are identical with the 
exception that the M&R method is much noisier. This verifies that this normal-state division is 
not necessary for S-Vac-S tunneling in this simple 'control' case of tunneling between similar 
elemental superconductors. This is an important statement for future experiments. Clearly it 
is experimentally simpler to take one measurement and use raw output rather than to take two 
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FIG. 3: (a) Comparison of occupied —d 2 I/dV 2 data by HHS experiment (blue), (J^Isc/^In) exper- 
imental data (black), Eliashberg calculation based on model a 2 F(co) (green), (b) Model a 2 F(co) altered 
from the Lorentzian model of Scalapino et al. ||29| . (c) Experimental dl/dV curve overlaid with Eliashberg- 
based simulation (green). Note the strong agreement with the exceptions of the large dips symmetric about 
£f (white arrows). 
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measurements at different times, under different conditions, and in combination with a numerical 
derivative combining the data sets. We have now shown that this efficient single HHS measurement 
is sufficient. (It is important to note that this is valid since the normal-state spectra is essentially 
flat. In the specific case of the hole-doped cuprates, where the normal state is highly non-trivial 
and spatially dependent, such a division might still be necessary |jT2~l| ). 

Since our spectra is so similar to that of M&R it should be straightforward to calculate the 
tunneling spectra with Eliashberg theory and match it to our experiment. The tunneling cur- 
rent between two superconductors separated by an insulating (or vacuum) barrier is given by 
I(V) oc Jq dcoNsc,i((o)Nsc,2((0 + V)[n F (co) -n F (co + V)] where V is the sample bias, n F (co) is 
the Fermi occupation number, and Nsc,i((°) an d ^sc,2(w) are the density of states (DOS) for the 
two respective superconductors. Once I(V) is known, the differential conductance dl/dV and the 
second derivative d 2 I/dV 2 can easily be found numerically. For a strong-coupling superconductor 
such as Pb, the DOS is given by N S c((o)/N N (0) = Re[w/ a/ co 2 - A 2 (ct>)] , where the S (N) sub- 
script denotes the superconducting (normal) state, and N^(0) is the DOS at £p. Here, A(co) is the 
complex frequency-dependent gap function given by the ratio of the anomalous self-energy and 
the renormalization function A(fi)) = <j)((o)/Z((o). Z(co) and (j>((o) are obtained by solving the 
real-axis Eliashberg equations Il29ll30ll using a 2 F{(o) as an input and our calculations are there- 
fore dependent on this choice. In this case we have two phonon spectra to consider, one for the 
bulk sample and one for the tip. For the former we follow Scalapino et al. [29] and define a 2 F 
with a sum of Lorentzians (for more details see [SM]). For the latter we modify this model with 
additional peaks to account for the thin-film nature and geometry of the superconducting tip where 
coupling can occur with phonon branches normally silent in bulk Pb. Figure [3] shows these model 
a 2 F functions and the resulting dl/dV and d 2 I/dV 2 which agree with the experiment (acquired 
in the two previously described experimental methods). The strongest difference between the two 
are the strong dips immediately past the coherence peaks in dl/dV [indicated by the arrows in 
Fig.^c)]. These dips are also manifested in d 2 I/dV 2 where the curves split at low energies. 

The pronounced dips in dl/dV are visually reminiscent of other tunneling studies — for exam- 
ple, break junction tunneling experiments on hole-doped BSCCO-2212 OTTl where the authors 
claim a connection between the spectral dip shapes (often referred to as peak-dip-hump) and 
bosonic modes. Similar peak-dip-hump structures have been seen numerous times in N-Vac-S 
standard STS studies on hole-doped BSCCO-2212 as well. Here we make no statement on the 
origin of these dips in the high-r c superconductors. We simply note that in our elemental super- 
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conductor study very similarly shaped dips occur and, in this benchmark case, their absence in 
a 2 F((Q) excludes those dips from being relevant to pairing. Similarly shaped dip features have 
previously been seen in STS S-Vac-S junctions Il25ll32ll33l , but none of these studies give a com- 
pelling explanation for their existence. None of those studies though were able to truly differenti- 
ate that dip feature from other bosonic signals as in this work. The most prevalent hypotheses for 
these features invoke proximity effect where extra features in d 2 I/dV 2 are observed in combined 
S-I-NS junctions lT34ll . surface impurity considerations lT32ll . or ^-selective tunneling processes in 
two-band superconductivity [|33ll . While a multiband model might not be appropriate for the case 
of Pb, the underlying idea of the ^-selective tunneling might be relevant due to Pb's anisotropic 
superconducting gap [|35l[36l . In this experiment the direction of anisotropy for the tip might not 
align with the anisotropic axis of the bulk. This alignment mismatch may relate to the similar 
dips in the break junction results between superconductors but does not elucidate peak-dip-hump 
behavior in standard STS studies with a normal metallic tip. 

The sandwich configuration of the M&R device necessarily demanded that occupied and un- 
occupied states be equivalent since the definition of the "sample" electrode and "tip" electrode 
could always be re-designated due to identical geometry of the two electrodes. In this experiment, 
the tunneling electrodes are distinguishable allowing us to differentiate between the occupied and 
unoccupied states of the bulk sample. Due to the very weak particle-hole asymmetry in Pb, theory 
demands — and our simulations demonstrate — that no differences between the occupied and unoc- 
cupied sides should exist. In our measured spectra the two sides are indeed mostly similar, but 
variations are in fact present. This particle-hole asymmetry does not have an obvious explanation. 
Immediately, we can discount the notion of this being an artifact of the tip structure. While it 
is possible for the electronic structure of the tip to add slope information or distort the measured 
spectrum through convolution, there is no present theory to account for the tip subtly modifying 
the phonon structure on only one side of £p. To explore these differences more thoroughly, both 
occupied and unoccupied d 2 I/dV 2 spectra were fit with Gaussian peaks. Two different fitting 
schemes were explored. As a starting point in fitting, the individual peak positions were con- 
strained to be identical for both the occupied and unoccupied states. The experimental asymmetry 
could be captured by either introducing an extra peak at Q. = 2. SI meV in the unoccupied side 
or by allowing the FWHM of the dual peak structure near 4 meV to have different values for the 
occupied and unoccupied sides (see data in Fig. 2). Neither of these fits have a physical basis 
within Eliashberg theory, however — extra phonon structure cannot exist on only one side of £p, 
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FIG. 4: Calculated a 2 F(a>) from spectral inversion of d 2 I/dV 2 treating the occupied (blue) and unoccupied 
(magenta) electronic states as separate spectra. The yellow curve is the device a 2 F(co) from M&R (H. 
Arrows mark features discussed in the text. 

nor is there a physical reason for the broadening of the phonon peaks to have different values on 
alternate sides of £p. Even more puzzling is the fact that only certain phonon modes seemed to be 
broadened asymmetrically about £p while the other modes have broadening symmetric about £p. 

We now advance from the forecasted model of a 2 F((o) to a full inversion of the experimental 
data to calculate a 2 F((d). Since the Eliashberg -based inversions are not compatible with occu- 
pied/unoccupied asymmetry calculations, this asymmetry can only be explored by treating the 
occupied and unoccupied sides of the spectra as independent data sets. Another complexity of this 
calculation is the question of how to treat a 2 F{(o) in terms of the sample superconductor and tip 
superconductor. Since the inversion process is based on an iterative process to bring the calculated 
d 2 I/dV 2 progressively closer to the experimental d 2 I/dV 2 , leaving both the tip and the sample as 
separate fitted spectra would result in an infinite number of different solutions to a 2 F{(0). The 
most impartial way then to perform the inversion is to constrain the magnitudes of a 2 F(ft)) t i p and 
a F(cq) sample m a fixed ratio determined by their respective A values while allowing their linked 
overall shape to vary. (For a thorough description of the fitting procedure see [SM]). 

Fig. |4] shows the inverted OC 2 F((o) occup i e( i and OC 2 F(co) unoccU p ied and its general agreement with 
the M&R result. Using the inverted a 2 F((o) spectrum we now have the information necessary to 
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extract the dimensionless coupling strength X = Jq 2a 2 F(o)) da/ a from our experimental data. 
The occupied and unoccupied curves give A occup ied = 1-38 and Aunoccupied = 1-34 with Average = 
1.36 compared to the M&R result of Ad ev ice = 1-33 ITJ. (From a particle-hole symmetric theory it 
is again unphysical to quote that the coupling strength and thus the T c is variable for the occupied 
and unoccupied sides of the spectrum, but from an experimental point of view we have two valid 
data sets to be used to obtain A values). 

The two strongest differences between this experiment's calculated a 2 F((d) and the M&R 
result are the additional spectral weight around 7 meV (Fig.|4} green arrow) and the splitting of the 
first peak (Fig.|4} white arrow). These can easily be explained by the thin-film nature of the Pb tip. 
In a first-principles phonon simulation [371 it was seen that in two-monolayer films the surface 
lifted the degeneracy of the transverse mode resulting in distinct transverse peaks. It was also 
seen that increased electron density on the surface resulted in stronger Pb-Pb bonds, stiffening the 
longitudinal peak which accounts for the added spectral weight seen around 7 meV. (The splitting 
of the transverse modes can also be explained by enhanced spin-orbit coupling [|38ll . but this option 
does not explain the weight at 7 meV). 

This then takes us to the differences between the occupied and unoccupied inversions of 
a 2 F{(0). The black arrows in Fig. [4] highlight these differences. The largest difference is an 
additional peak in the transverse mode (black arrow near 5 meV). This peak sits roughly in the 
middle of the range of 3-7 meV in which the occupied and unoccupied spectral data display the 
most significant contrast (Fig. 2), and is near the 5-meV energy where some shear vertical surface 
phonons have been measured at energies that have been potentially softened by spin-orbit cou- 
pling 113911 . This does not offer any explanation as to why this would cause this phonon signal to 
be different when measured on the occupied vs. unoccupied side. The other difference between 
the occupied and unoccupied inversions is that a F(c&) occupied has two very small peaks, one near 
2 meV and one near 10 meV. Since the coupling strength scales inversely with ft), the 2-meV peak 
is of far greater importance. A peak near this 2 meV energy is also seen both in the occupied and 
unoccupied spectra of d 2 I/dV 2 (although at different strengths). This peak energy is very close 
to the Kohn anomaly in the transverse acoustic phonon at q = (0.26,0.26,0) measured at energy 
E = 2.32 meV recently by inelastic neutron scattering [40] . It had been suggested that such a 
Kohn anomaly is caused by spin-orbit coupling or many-body interactions and has a correlation to 
the size of the superconducting gap 2A. This correlation though is controversial and may be coin- 
cidental flu . Nothing about the way we understand these anomalies lends a suitable explanation 
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as to why they would appear stronger or weaker (or not exist at all) in the occupied vs. unoccupied 
spectra. Theories that include electron-hole asymmetries [42] may need to be considered. 
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